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ABSTRACT 

The radio signature of a shock travelling through the solar corona is known as a type II solar radio burst. In rare cases these bursts 
can exhibit a fine structure known as ‘herringbones’, which are a direct indicator of particle acceleration occurring at the shock 
front. However, few studies have been performed on herringbones and the details of the underlying particle acceleration processes 
are unknown. Here, we use an image processing technique known as the Hough transform to statistically analyse the herringbone 
fine structure in a radio burst at ~20-90MHz observed from the Rosse Solar-Terrestrial Observatory on 2011 September 22. We 
identify 188 individual bursts which are signatures of bi-directional electron beams continuously accelerated to speeds of 0.16^^ }q c. 
This occurs at a shock acceleration site initially at a constant altitude of ~0.6 R© in the corona, followed by a shift to ~0.5 R©. The 
anti-sunward beams travel a distance of IVO^^^^^Mm (and possibly further) away from the acceleration site, while those travelling 
toward the sun come to a stop sooner, reaching a smaller distance of 112^^^ Mm. We show that the stopping distance for the sunward 
beams may depend on the total number density and the velocity of the beam. Our study concludes that a detailed statistical analysis 
of herringbone fine structure can provide information on the physical properties of the corona which lead to these relatively rare radio 
bursts. 
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1. Introduction 


The longest known signature of shocks propagating into the so- 
lar corona are t ype II radio bursts ( |Wild et ^ |1959[ [Nelson 
[Melrosej |1985| ). It has been postulated that they are produced 
by a disturbance travelling into the solar atmosphere at speeds 
of over 500 km s“^ About 20% of type IIs contain a fi ne struc- particles, 
ture of pulsations or ‘thorns’, known as herringbones ([Roberts | 

1959[), wh ich typically occur between 10-120 MHz ( [Cairns 


(it is generally the opposite in type Ills), and the sense of polar¬ 
isation of their fundamental and harmonic bands is opposite to 
that of type Ills. Indeed, later findings eventually confirmed a 
closer link to shock activity in the corona, with more intense 
type IIs being more likely to have herringbone structure (Cane 
& White|[1989[). Regarding the pr oper ties of shock accelerated 


Mann & Classen 


( [1995[ ) and [Mann & Klasseli] ( |2005[ ) 


Robinson| 1987[ ), and are thought to be observations of electron 


beams escaping the shock as it propagates. Although herring¬ 
bone radio bursts provide a wealth of information on the particle 
acceleration process inside coronal shocks, few studies of their 
properties have been made. A statistical analysis of herringbone 
features could yield important information on the type of accel¬ 
eration process occurring in coronal shocks and/or the environ¬ 
ment into which particles are accelerated. 

Observation s of herringbone bursts were first report ed by 
Roberts[ ([1959[ ) and a description was later provided by Wild[ 
et al.[ ( [1959| ) whereby 'shock fronts propagate normal to the 


showed that the electron beams that cause herringbones have 
speeds of ~0.1 c. [Carley et al.[ ( [2013| ) reported similar speeds 
and showed that the acceleration process has a quasi-periodicity 
of between 2-11 seconds and may be due to particle acceleration 
at an expanding CME flank. 


magnetic field and continuously eject bunches of fast electrons 
along the field '\ this assertion was also illustrated in [Stewart &[ 
[Magun[ ( [1980[ ). Although herringbones appear quite similar to 
bi-directiona l type Ills, such as those rep orted in [Aschwanden[ 
[et al.[ ( |1995[ ), [Cairns & Robinson] ( [1987[ ) concluded that her¬ 
ringbones are a different phenomenon to type III bursts; their 
frequency-time profiles are characteristically different, -80% of 
herringbones have stronger fundamental emission than harmonic 


There has been little theoretical investigation as to the origin 
of the bursts. Eor example, what gives herringbone radio bursts 
this quasi-periodicity or ‘burs tines s’ in tim e? And what are the 
coronal s tructures that lead to herrin gbones ? [Zlobec et al.[ ( fl993[ ) 
and later [Vandas & Karlicky ( |2011[ ) suggested that particles ac¬ 
celerated in a magnetic trap on a wavy shock front might be 
responsible. Similarly, it has been suggested that the burstiness 
is due to inh omogeneity and repeated electron acceler a tion on a 
wavy shock (|Guo & Giacalone| [2010[ [Burgessj [2006| ). [Schmidt 
& C airns I ( [2012| ) presented such a scenario in a model of a rip¬ 
pled shock on a CME flank accelerating particle beams into 
the surrounding corona. Alternatively, herringbones may be pro¬ 
duced in the termination shocks of super-Alfvenic outflow jets of 
a reconnection region ( [Aurass et al. [2002[ [Mann et aT} [2009[ ). 
Such a theory may explain why some herringbones appear to 
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come from a stationary height in the corona (Aurass & Mann 
[2004] ). 

Despite the combination of observation, theory, and mod¬ 
elling that has gone before, there is no complete explanation 
of herringbone fine structure. An understanding of these bursts 
could benefit from a st atistical analysis such as was previously 
perfor med for type Ills ( Saint-Hilaire et al.[|2Q13t[Lobzin et al 


|2Q11| ). However, because of the high number of individual her¬ 
ringbones in a typical radio burst, an automatic feature recog¬ 
nition technique may be required to stud y them accurately. One 
such technique is the Hough transform ( Hough[ 19611 |Duda & 
Hart[|1972] ), used to identify straight lines in images. This tool 



was employed as a feature recognition te chnique for type Ills 
as part of the HELI O Feature Catalogue ( [Bentley et al.[ [201 1[ 
Bonnin et al.||2011|) and also for automatic recognition of type 


10:00 10:20 ^1^:40 11:00'-^ 11:20 11:40 12:00 

^ ^ Start Time (2011-Sep-22 10:00^00OT)^ 


RSTO CALLISTO ' 


IIs ( jLobzin et al.||2010| ). Those studies were for large data sets; 
here we constrain ourselves to a single radio event exhibiting 
herringbone bursts to test the applicability of the Hough trans¬ 
form to this kind of radio fine structure. 

Section 2 concentrates on radio observations of the event; 
Section 3 provides a description of the Hough transform and 
methods; Section 4 gives a statistical analysis of herringbone 
burst properties; and Sections 5 and 6 provide a discussion of 
the physical interpretation of the bursts and conclusions, respec¬ 
tively. 

2. Observations 

The SOL201 l-09-22T10:30 event (XI.4, Figure[^) was accom¬ 
panied by a variety of radio activity including type IIs, type 
Ills and herringbone radio bursts. The fiare occurred at N09F89 
and was accompanied by a CME, EUV wave and a number of 
radio sources imaged by the Nan 9 ay Radioheliograph (NRH; 



10:51:00 10:51:30 10:52:00 10:52:30 10:53:(j0 

Start Time (2011 -Sep-22 10:50:30 UT) \ 

Back Projected Hough Transform ^ 



Kerdraon & Delouisl 

1997 

); a detailed analysis of the event is 

10:50:30 10:51:00 10:51:30 10:52:00 10:52:30 10:53:00 

given in [Carley et al.) 

2013 

). The herringbone radio activity be- 

Start Time (2011-Sep-22 10:50:30 UT) 


gan at -10:47:30 UT between 10-90 MHz, detected using the 
lowest frequency Cal listo receiver at the Rosse Solar-Terrestrial 
Observatory (RSTQ; [Zucca et al.]|2012| ), at Birr Castle, Ireland. 
Herringbones drifting to both high (reverse drift) and low fre¬ 
quency (forward drift) can be seen in Figure [^. An interesting 
feature of this burst was the lack of drift of the burst ‘back¬ 
bone’ i.e. the herringbones appear to originate from t he same 
frequency over time, similar to observations reported in [Aurass 
|er^ ( |2n5^ ; [Mel’nik et al] ( [20041 ); [Mel ’Nik et al.[ ( |2005[ ). This 
is suggestive of particle acceleration from a shock at a constant 
altitude, or at least in an environment of constant density in 
the frame of the shock. Usually, the backbone drifts to low fre¬ 
quencies like a normal type II while emitting herringbones (see 
Miteva & Mann ( [2007 [ ) for a good example). 


A previous analysis showed that there is some level of peri¬ 
odicity of be tween 2-11 seconds on the occurrence of the bursts 
in this event ( [Carley et~SH 2013[ ). Our goal here is to further sta¬ 
tistically analyse the bursts in terms of drift rate and intensity as 
well as the velocity and distance travelled of the electron beams 
causing the radio emission. 

3. Method 

To perform a statistical analysis we used the Hough transform 
( [Hough[ [1961[ ), a feature recognition algorithm that is used to 
identify straight lines in images. In the Hough transform any 
point in an image may be represented in a polar coordinate 
‘Hough space’ as a sinusoid; if two points in the original im¬ 
age lie on a straight line, their sinusoids intersect in the Hough 


Fig. 1. (a) GOES soft X-ray light curves showing an XI.4 class 
fiare on 22 September 2011 at ~ 10:30 UT. (b) During the im- 
pulsive/peak phase of the fiare the RSTO 10-100 MHz Callisto 
receiver observed herringbone radio bursts from -10:48- 
10:55 UT; only a section of all herringbone activity is shown here 
for clarity, (c) The back projected Hough transform of the area 
containing the reverse drift bursts. This data contains positions 
of all bursts, which are much more clearly defined and subject to 
less RFI or background noise. We note that the Hough transform 
provides the positions of burst maxima but does not produce the 
stop frequency in the dynamic spectra. For our analysis the stop¬ 
ping frequency was chosen manually. For a complete overview 
of all radio activity from this event, see [Carley et al.[ ( [MT3 ). 


space. Lines in the original image space are found by searching 
for the intersec tion points of sinusoids in Hough space, see [Duda[ 
& Hart \\912) for more details. If the image is transformed into 
Hough space, it is possible to ‘back-project’ the Hough space to 
re-obtain the original image. Furthermore, it is possible to back- 
project a particular region of the Hough space, such that only 
lines with a particular orientation in the original image will pass 
through the transform. This is very useful for filtering an image 
for linear features of a particular orientation or slope. 

For the herringbone observations, we firstly pass our dy¬ 
namic spectrum through the IDL gradient.pro function to pro¬ 
duce an image (dynamic spectrum) of intensity gradient. This 
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eCallisto, Birr Castle, Ireland. 



10:51:40 10:52:00 10:52:20 10:52:40 10:53:00 

Start time: 2011-09-22710:51:29 (UT) 


Fig. 2. (a) Indication of each burst detected in the algorithm 
using the Hough transform. The forward drift bursts are high¬ 
lighted in orange-red and reverse drift are in green-blue. In gen¬ 
eral, the forward drift were more difficult to detect due to their 
weaker signal to noise and because they drift into an area of 
high RFI <20 MHz. This is the reason why there are more re¬ 
verse drift detections than forward drift, (b) Zoom of the second 
set of reverse drifting bursts with detections overlaid as circles. 



Fig. 3. Electron density map of th e solar corona prod uced from 
AIA and LASCO C2 observations ( |Zucca et ^|2014| ). A density 
profile of the corona for the data of the event was produced by 
averaging the density as a function of radius between the lines at 
position angles of 90° and 135°. This density as a function of ra¬ 
dius was used in deriving electron beam kinematics from the her¬ 
ringbone radio burst statistics. The inner and outer black points 
represent the heights at which 43 MHz (0.47 R©) and 32 MHz 
(0.59 R©) are expected to occur. 


is to produce an edge detection of the radio bursts, delineating 
the burst peaks more sharply. Next, we choose a section of the 
dynamic spectrum containing only the reverse drift bursts, such 
as the dashed box in Figure [^, and pass it through the Hough 
transform. Since all bursts in this section have a particular ori¬ 
entation, we know generally which area of Hough space from 
which to produce a back-projection. This allows us to re-obtain 
the original image but only containing these bursts specifically, 
see Figure [^. Having been passed through the Hough transform 
and back-projected, the bursts are now sharp, unbroken lines 
which are much more clearly defined and subject to less RFI or 
background noise. These sharply defined peak positions allow 
us to identify the burst peak for each frequency in the dynamic 
spectrum using a peak finding algorithm. 

We note that the Hough transform fails to demarcate the start 
and end frequencies of the bursts. In the back-projection, the lin¬ 
ear features can continue after the burst has stopped. For exam¬ 
ple, by comparing Figure and Ic we can see some herring¬ 
bones drift no higher than 60 MHz, but in the the back-projection 
all detections extend as far as 90 MHz. This is because the trans¬ 
form works by assigning a line to two or more points in the orig¬ 
inal image (dynamic spectrum), this line is represented mathe¬ 
matically in the transform without bounds. In order for a burst 
end frequency to be chosen, a user must go through the bursts 
one by one and select an end frequency i.e. where the burst fi¬ 
nally fades into the background emission. This ensures that the 
correct burst length in frequency is chosen. 

The total amount of detections is shown in Figure with 
forward drift bursts in orange-red and reverse drift in green- 
blue, with bottom panel showing a zoom for detail. The detec¬ 


tion performs well, detecting 188 bursts in total. However there 
are a number of shortcomings. There is not 100% accuracy, and 
the transform fails to detect bursts in a particularly noisy back¬ 
ground environment, or in regions of low signal to noise. This 
is particularly notable in a region midway through the activity at 
10:50:30 UT during which there is a shift of the backbone from 
~32 MHz to ~43 MHz (we discuss the physical interpretation of 
this in Section 4). During this shift, the activity becomes quite 
‘patchy’ and bursts are not so easily detectable, especially the 
reverse bursts. This is why detections are sporadic around this 
area. Also, the forward drift bursts (those drifting to lower fre¬ 
quencies) are much weaker and generally drift into an environ¬ 
ment of high RFI in the AM band below ~20 MHz, hence de¬ 
tections of forward drift bursts were less successful than reverse 
drift. In fact, after the backbone shift, the forward drift bursts 
became very weak, and not easily detectable. These few bursts 
(yellow detections. Figure [^) had to be selected by point-and- 
click. A test of detection threshold in a quiet area of our dynamic 
spectrum revealed that detection is successful if the burst has an 
intensity greater than ~3(T above the quiet background, where 
cr is the standard deviation of values in the quiet background. If 
the burst is surround by noise of comparable intensity (RFI, for 
example), detection is much more difficult. 

There is a region around 10:48:00-10:49:30 at 60-80MHz 
(dashed box. Figure [^) where features were left unselected. 
They are much fainter and more diffuse than the herringbones, 
with no discernible drift. We consider these to be possibly part 
of a continuum emission, so they are excluded here. 


3 

































E. Carley et al.: Radio Observations of Bi-directional Electron Beams 



Fig. 4. (a) Frequency vs time data points for all detected 
bursts for both reverse (blue-green) and forward drifting bursts 
(orange-red). One line represents one burst, with 188 in total, (b) 
Linear fits to each burst detection. The green-red lines originate 
from ~32 MHz and are detections of the first set of bursts be¬ 
tween 10:47:30-10:50:30 UT. The blue-orange lines are the sec¬ 
ond set of bursts originating from ~43 MHz between 10:50:20- 
10:53:00 UT. There are not equal numbers of positive and re¬ 
verse drift bursts, indicating there is not a one-to-one correspon¬ 
dence of particle beams travelling toward and away from the sun, 
i.e. beams are bi-directional but not always in reverse and for¬ 
ward drift partners. 


Overall, the algorithm is successful in identifying individ¬ 
ual herringbone structure for the majority of time and frequency 
ranges. It is essentially a semi-automated routine, whereby the 
transform detects all possible bursts, with the user providing 
some input to evaluate the validity of the detection. 


4. Results 

In the following we refer to the herringbones as observed in dy¬ 
namic spectra as bursts, while the population of fast electrons 
causing the radio burst are called a ‘beam’ i.e. one beam of elec¬ 
trons causes one herringbone burst. 

The detection of 188 bursts allows for the statistical analysis 
of some useful properties. The total number of reverse and for¬ 
ward drift bursts are 118 and 70, respectively. It is interesting to 
note that there is not a reverse drift for every forward drift burst 


and vice versa. This shows that the beams are bi-directional in 
space, but do not always emerge from the shock as sunward and 
anti-sunward partners (as noted in |Caims & Robins^ ( |1987| )) - 
we discuss this further in Section 5. 


In the following we derive kinematical results of the particle 
beams causing the herringbones. To do this we make the assump¬ 
tion that we observe fundamental plasma emission. It is reason¬ 
able to assume that we observe the fundamental component of 
the burst, given that fundamental fiux exceeds harmonic fiux in 
80% of herringbone observations ( |Caims & Robinson 1987 ). 
However, as we discuss below the general results for beam kine¬ 
matics do not change if we assume harmonic emission. 


4.1. Burst altitude and backbone shift 


It is apparent from Figure that mid-way through the radio ac¬ 
tivity at -10:50:30 UT, there is a shift of the backbone to higher 
frequency, from -32 MHz to -43 MHz. This could indicate the 
acceleration region shifted in altitude midway through the radio 
activity. To calculate this altitude shift we firstly assume funda¬ 
mental plasma emission (fplasma) and convert to electron number 
density (Ue) via fplasma ~ 8980 The height at whi ch this 


density occ urs was found using a TCD density model ( |Zucca 
|et al.||2014] ), produced from density measurements of the corona 


using the Atmospheric Imaging Assembly (AIA; [Lemen et al. 


2012|) and Large Ang le Spectrometric Coronagraph (LASCO; 
Brueckner et al.[p^95]) C2 observations on the date of the event. 


see Figure |3| An average of the density profile was taken be¬ 
tween position angles of 90° and 135°, i.e. where this event took 
place off the limb (see |Carley et ar] ( |2013| )). This density profile 
(produced from observations) prevents the needs for a ‘guess’ 
model, and allows us to be more confident that the density model 
used in the analysis is a correct description of the corona for this 
event. We find that the acceleration region is initially at constant 
altitude of 0.59/?© then undergoes a significant shift to a lower 
altitude of 0.47 Rq. 


Although a shift in frequency is generally taken as an altitude 
shift, we cannot rule out that a shift in frequency may also be 
from an accelerator at a constant altitude encountering a change 
in background density profile. For example, if the accelerator 
propagated at constant altitude, but changed in position angle 
it may encounter larger densities, resulting in a shift to higher 
emission frequencies for the herringbones. Figure shows that 
different density profiles exist within the bounds of where the 
event took place. We stress that the shift in frequency is due to an 
accelerator moving in an inhomogenous corona, either shifting 
in height or position angle (or both). Either way, it is clear that 
the accelerator changes environment midway through the burst. 
We show later in Section 4.3 that this results in a change in burst 
intensity profile and speed. 

It is also worth noting the density calculated from the TCD 
model is produced from a line-of-sight average. This inherently 
ignores the spatial variability in the corona, for example at the 
interface between open and closed field. The actual densities 
may be h igher or lower than the average we have calculated. 
However, |Zucca et ak] \2b\4) showed that for this event the ex¬ 
pected heights we calculate from the density map are a close 
match to the heights of radio sources observed by NRH, con¬ 
firming that our density values are reliable. 


4 



























E. Carley et al.: Radio Observations of Bi-directional Electron Beams 



-10 -5 0 5 10 15 

Drift rate (MHz s') 



0.0 0.1 0.2 0.3 0.4 0.5 

Beam displacement (Rq) 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 

Velocity (c) 


Fig. 5. (a) Histogram of drift rates of reverse drift (blue stripes) 
and forward drift (red) bursts. The reverse drift bursts have a 
bimodal distribution. Those with a generally higher drift rate 
(above lOMHzs"^) are the reverse drifting bursts that start at 
43 MHz. Taking the absolute value of these drifts and using a 
TCD density model we produce a histogram of velocities for the 
electron beams causing the herringbones (b), with an mean of 
0.16 c. We note that although the drift histogram is bimodal, the 
velocity histogram is not. This is due to electron beams of a sin¬ 
gle velocity distribution propagating in different heights of the 
atmosphere. A different height means a different density gradi¬ 
ent and hence a different drift rate. 


4.2. Burst drift rates and particle beam speeds 

Figure shows the frequency time traces of all detected bursts 
(a), along with linear fits to the data (b). These linear fits were 
used to obtain the drift rates of all herringbones in the radio burst 
and produce a histogram, as shown in Figure [^a). The bin size 
is calculated from the smallest shift detectable in the dynamic 
spectrum, taken to be two pixels in frequency (0.9 MHz) and 
two pixels in time (0.5 seconds) giving a minimum detectable 
drift of 1.8 MHz s“^ The histogram of burst drift rates were 
also converted to beam velocities using the TCD density model 
Figure [^b). The bin size for velocity is the minimum detectable 
drift rate converted to velocity using our density model (0.04 c). 

The drift rate histogram has some notable features. Forward 
(red) and reverse (blue stripes) drifts have similar absolute 


Fig. 6. Histogram of beam displacements, produced by convert¬ 
ing the frequency span of individual bursts in the dynamic spec¬ 
tra to distances using a TCD density model. An interesting re¬ 
sult here is that the forward drift bursts (those travelling away 
from the sun) tend to travel larger distances than the reverse drift 
(those travelling toward the sun). It is clear that the positive drift 
bursts solely make up the high-end tail of this histogram. This 
may be an indication that electron beams travelling toward the 
surface stop travelling (or stop emitting) sooner than those trav¬ 
elling away from the sun. The bin size is calculated from mini¬ 
mum detectable shift over two pixels (0.9 MHz) in the dynamic 
spectrum, giving a 0.03 Rq bin size. 


magnitudes, however the reverse drift bursts appear to be bi¬ 
modal. The high drift population (greater than ~10MHzs“^) 
belong to the reverse drift bursts which start at 43 MHz (be¬ 
tween 10:50:20-10:53:00UT). This is not necessarily because 
the beams causing these bursts have a higher velocity, but be¬ 
cause they are accelerated lower in the atmosphere, and therefore 
encounter a steeper gradient in density, giving them a steeper 
drift in frequency space. This can be confirmed by analysing the 
velocity histogram in Figure [^b). Forward and reverse beams 
belong to similar and overlapping velocity distributions, with 
their total histogram peaked at ~0.16c (maximum and mini¬ 
mum speeds of 0.36 and 0.06 c, respectively , similar to those re- 
ported for herringbones at -10-120 MHz by |Caims & Robinson 
( 1987| )), with the beams producing forward and reverse drift 
bursts showing only a slight difference in speed distribution. This 
is a good indication that we have detected drift rates accurately, 
i.e. independent drift rates and varying start frequencies result 
in similar velocities for all electron beams. This is evidence that 
both forward and reverse beams throughout the entire radio burst 
belong to the same accelerator even though there is an altitude 
shift halfway through the activity. The effect of assuming har¬ 
monic emission does not change this result, it simply shifts the 
velocity histogram to higher values by -0.05 c. 


4.3. Beam displacements 

Figure ^shows a histogram of electron beam displacements (dis¬ 
tance that the electron beam travels) in space. The displacements 
were calculated by converting start and end frequency to start 
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Table 1. Compilation of radio burst properties and derived particle beam kinematics using a TCD density model. F is forward drift 
and away from the sun, while R is reverse drift and toward the sun. 



# 

Drift (MHz s-^ 
min mean max 

Time (s) 
mean mode 

Vel. (c) 

min mean max 

Displ. (Mm) 
min mean max 

Stop altitude (R©) 
min mean max 

R 

118 

3.0 8.1 14 

2.8 2.1 

0.06 0.14 0.25 

36 112 197 

0.2 0.34 0.45 

F 

70 

-8.5 -4.6 -2.8 

3.3 2.6 

0.11 0.18 0.27 

74 170 344 

- 


and end height using our density model. The bin size is calcu¬ 
lated from the minimum detectable frequency shift in the dy¬ 
namics spectrum (2 pixels = 0.9 MHz) converted to a displace¬ 
ment from the model (0.03 Rq). The overall distribution appears 
bimodal, with the upward electron beams making up the high- 
end tail of the distribution (greater than 0.3 Rq). 

This is an interesting result, showing that the beams trav¬ 
elling toward the solar surface are either stopped sooner or 
stop emitting sooner. Sunward beams travel in the range of 36- 
196 Mm with an average of 112 Mm, while those travelling away 
from the sun have a larger displacement in the range of 73- 
344 Mm with an average of 170 Mm. In fact, the upward-beam 
displacements may be longer than this; the stopping frequency of 
forward drifting bursts is more difficult to identify, as the bursts 
generally drift into a region of high RFI and below ~20 MHz 
toward the ionospheric cutoff frequency. This is the reason for 
the bi-modality of the forward drift bursts. Those forward drift 
bursts that belong to the small displacements in the distribution 
(-0.15 Rq in the red histogram) primarily belong to the first set 
of forward drifters in Figure [^. It can be seen here that some of 
these detections (red circles) are cut short because they drift into 
a region of RFI, just below 30 MHz. These bursts likely drift to 
lower frequencies and hence should have larger displacements. 
The reason for upward beams travelling further (or emitting for 
longer) than sunward beams is discussed in the Section 5. 

Since the velocity of both reverse and forward drift parti¬ 
cles are the same, the lifetime of the forward particles should 
be slightly longer, given that their displacements are longer. We 
confirm this by looking at the mean and mode lifetime of the 
bursts in Table 1, which shows that the forward particles gen¬ 
erally have a longer lifetime by -0.5 second, i.e. travelling for 
0.5 extra second at speeds up to -0.27 c gives an extra distance 
of -40 Mm, which is similar to the extra displacement range we 
find for the forward bursts. 

Again, the effect of assuming harmonic emission does not 
change the result. Both forward and reverse drift beams are 
shifted to smaller displacements by -15 Mm, however the for¬ 
ward drift beams still make up the high-end tail of the displace¬ 
ment histogram and on average travel a further distance. The 
radio burst and derived electron beam kinematic properties, as 
well as the difference between forward and reverse drifters, are 
listed in Table 1. 


Finally, we tested all kinematics results using the popular 
Saito, Leblanc, Newkirk, Mann, and Baumbach-Allen density 
models (|Saito et aL| 1977t|Leblanc et ar| | 1998[[Newkirk||l96H 


Mann et al. 1999t Allen |1947 ), see Table Although the ab¬ 
solute values of the derived kinematics differ slightly, the gen¬ 
eral outcome is the same for all models. Table [2] also serves to 
highlight the variable results one might expect when using these 
density models to derive radio burst kinematics. Here we take 
the TCD model to be the most reliable, as it is based on density 
estimates from the day of the event. We reiterate that the TCD 
density model provides a reliable description of density values in 
the corona, given the close agreeme nt of NRH observe d source 
heights and the TCD model heights ( [Zucca et 5^12014| ). 


4.4. Burst intensity correlations 


In the following we analyse the correlations amongst burst in¬ 
tensity and other parameters, as well as the bursts rate of change 
of intensity over the burst lifetime. Figure [TJa) gives a single 
example of how a burst maximum intensity changes over its 
lifetime, or equivalently, changes with respect to frequency. In 
Figure iTjb), we plot burst drift rate vs the rate of change of in¬ 
tensity over the burst lifetime. This plot contains only the reverse 
drift bursts. This was done because the reverse drift bursts have 
a much better signal to noise, are subject to much less RFI and 
their intensity can generally be sampled for most of their dura¬ 
tion - this is not the case for the forward drift bursts. 

The system at RSTO is uncalibrated, so intensity is quoted 
in data numbers (DN). We can only note that the receiver re¬ 
sponse is logarithmic in intensity. The data have also been back¬ 
ground subtracted using the SolarSoft constbacksub.pro routine. 
This ffnds a relatively quiet spectra of the spectrogram (contain¬ 
ing no bursts) and subtracts the quiet spectra from each time 
step in the spectrogram. Since the RSTO Callisto sy stem has a 
logarithmic response to any signal ( |Benz et aT 2005| ), a subtrac¬ 
tion of a background in this manner is essentially a division of 
the spectrogram by the background. This flattens any variation 
in intensity that the system may have across the bandwidth, so 
changes in intensity due to instrumental response may be ruled 
out. Furthermore, lack of knowledge of the absolute ffux values 
should not affect our general results. 

Figure [TJb) shows an overall negative trend with a correla¬ 
tion coefficient of -0.73. The negative trend reffects the fact that 
the bursts decrease in intensity over time, with the slope of the 
fit to the data being equal to dljdf i.e, the rate of change of in¬ 
tensity over the bursts’ duration in frequency. Decreasing inten¬ 
sity as the burst drifts in the dynam ic spectrum is the opposite 
to what is found in type III bursts ( [Cairns & Robinson| |1987[ 


Table 2. Comparison of beam kinematics for different density 
models. All quoted values are mean values. The beam displace¬ 
ments are shown for forward (F) and reverse (R) drifters sepa¬ 
rately. As expected, the values differ among the different models. 
However, the general results are unaffected by a change of den¬ 
sity model. 

t Only valid if the shift in frequency of the backbone is interpreted as an actual 
source altitude shift. 

* Baumbach-Allen. 


Model 

Altitude Shift t (R©) 

Beam Displ. (Mm) 

Beam Speed (c) 


start —> end 

F 

R 

F 

R 

TCD 

0.59 ^ 0.47 

170 

112 

0.18 

0.14 

Saito 

0.53 ^ 0.39 

218 

132 

0.22 

0.17 

Leblanc 

0.38 ^ 0.25 

202 

123 

0.2 

0.16 

Newirk 

0.74 ^ 0.57 

286 

141 

0.28 

0.18 

Mann 

0.35 ^ 0.25 

165 

90 

0.26 

0.17 

B-A* 

0.56 ^ 0.41 

261 

128 

0.16 

0.12 
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2 4 6 

Seconds after 10:51:38 UT 



Drift rate (MHz s^) 

Fig. 7. (a) Example of a single herringbone burst intensity at 
each frequency. The burst is most intense at its start and de¬ 
creases in intensity as it drifts to higher frequencies, i.e. de¬ 
creases in intensity over its lifetime, (b) Scatter plot of burst drift 
rate vs the rate of decrease in intensity along the burst for the re¬ 
verse drift bursts only. The receiving system is uncalibrated, so 
the intensity are simply quoted in Data Numbers or (DNs). The 
green points are the reverse drift bursts starting at 32 MHz (R 32 ), 
while the blue points are reverse drift bursts starting at 43 MHz 
(/? 43 ). We note that the R 32 have a lower drift and rate of change 
of intensity than the 7^43 bursts. Overall, a clear negative correla¬ 
tion is shown here; this occurs because the majority of the bursts 
decrease in intensity over their lifetime, i.e. the slope of the black 
line is {dlldt) = m{df /dt) => m = dl/dt x dt/df = dl/df. This 
is the opposite of type III bursts, which generally show an in¬ 
crease in intensity over their lifetime. 


[Saint-Hilaire et al.) |2013| ), meaning their may be a fundamen¬ 
tal difference to the propagation characteristics of type Ills and 
herringbones, or the environment in which the electron beams 
propagate may be different. 

Figure [7] also shows the separation of the distribution into 
two populations, those belonging to bursts that start at 32 MHz 
(green circles) and those that start at 43 MHz (blue circles). The 
green population has both smaller drift and smaller change in 
intensity, with some having an increase in intensity over time. 


Four of the bursts from the green population show a slow rise 
and fall of intensity with time. 

The different morphological characteristics of the 32 MHz 
bursts may be due to the altitude and/or coronal environment 
shift of the accelerator as described above. Once the acceleration 
region shifts to a new environment, the radio bursts increase both 
in drift rate and the rate of change of intensity (as seen from the 
change from green to blue points in Figure [7]). This suggests 
that the appearance of herringbones may depend on the kind of 
coronal environment in which the electron beams propagate. 

Finally, we investigate any correlation between intensity and 
beam displacement and velocity for the reverse drift bursts, 
shown in Figure [^a). The top figure shows beam displacement 
in space vs the maximum intensity of the bursts above the lo¬ 
cal background. The overall correlation is positive but weak, 
with a correlation coefficient of 0.48. However, the two popula¬ 
tions of reverse drift bursts again show a difference. Those with 
start frequency of 43 MHz (blue) are generally more intense and 
travel further than those starting at 32 MHz (green). The sepa¬ 
rate populations by themselves exhibit positive correlations ( 0.6 
and 0.55), showing the more intense the radio burst, the further 
the beam travels in space (toward the solar surface). We also test 
for any correlation between burst intensity and beam velocity 
shown in Figure [^b); we find almost no correlation here with 
correlation coefficient of 0.31. However we do note that the two 
populations of reverse drift bursts are separated, with the faster 
bursts generally being more intense on average. Again, the as¬ 
sumption of harmonic emission does not affect the above result. 


5. Discussion 

In this section we discuss the results in terms of the coronal 
environment of the herringbones and with regard to the stan¬ 
dard theories of radio burst intensity, i.e. wave mode conversion 
from Fangmuir waves to ion acoustic and electromagnetic waves 


(|Robinson| 

I992[|Robinson et al.||1994l|Fi et al.||2012l|Ratclilfe| 

let al.l|2014 



5.1. Burst starting frequencies 

An interesting feature of these herringbones is that both for¬ 
ward and reverse drift bursts originate at the same frequency 
(Figure [^a)). This may give clues as to the shock region in 
which the beam propagates and/or the beam formation and emis¬ 
sion mechanisms at play. 

Firstly, the reverse and forward drifters having the same start 
frequency indicates that electron beams travelling toward and 
away from the sun begin propagating in a similar density envi¬ 
ronment. This may be evidence that the forward and reverse drift 
bursts are not from beams ahead and behind the shock, respec¬ 
tively. If this were so, the reverse drift bursts would start at a sep¬ 
arate and higher frequency because of reverse electron beams en¬ 
countering a sudden downstream density jump across the shock. 
Instead we propose that forward and reverse drift bursts are both 
from electron beams accelerated upstream of the shock. For ex¬ 
ample, this would occur if a wavy shock propagated laterally in 
the corona and constantly accelerated electron beams onto field 
li nes in its upstream r egion, such as that illu strated in Figu re 5 
of |Zlobec et al.| ( p^93] ) or in the illustration of |Cliver| ( |2013| ). 

Secondly, we may interpret the similar start frequencies by 
comparing herringbones to type III observation and theory. The 
starting frequency of type III bursts depends upon a number of 
initial electron beam characteristics, namely the injection height. 
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size, time and energy distribution of the beam (Kane et al. 
IReid et al.| [MTT] |Reid & Ratcliffel |Reid et all 




The beam has to travel a certain distance before velocity dis¬ 
persion will create an unstable electron distribution, assuming 
the injection profile is a power-law. We do not observe any sep¬ 
aration in the starting frequencies between the forward and re¬ 
verse herringbone bursts. A lack of separation can be explained 
with very fast (< 0.1 s) injection times in an accelerator that 
is small along the dimension of electron propagation. This sce¬ 
nario could be feasible for electron acceleration originating in a 
small diffusive region in a shock. Another explanation is that the 
injected electron distribution is initially unstable to Langmuir 
waves. Essentially the accelerator creates a pre-formed beam at 
some velocity v » (where is the thermal electron tem¬ 
perature), leading to bi-directional bursts with no frequency gap. 
Such characteristics may be an indicator of the acceleration pro¬ 
cess at play in the shock. 

We stress that it may be upstream beam propagation andjor 
immediate beam formation that may lead to the same start fre¬ 
quency of reverse and forward drifters. 



Beam displacement (Rq) 


5.2. Displacements and emission mechanisms 


Although they have similar velocities, the forward and reverse 
drift bursts show a significant difference in the distances the elec¬ 
tron beams travel. The anti-sunward beams travel a distance of 
Mm away from the acceleration site (and possibly fur¬ 
ther, given the herringbones may reach frequencies below our 
observational limit at ~ 15 MHz), while those travelling toward 
the sun come to a stop sooner, reaching a smaller distance of 
Mm. The sunward beams either stop travelling and/or 
stop emitting possibly due to them encountering a different envi¬ 
ronment the lower in the corona they reach. The reason for this 
may be in the emission/absorption mechanisms at play in the 
corona, which has been investigated for type III bursts. 

Observationally, reverse type III bursts are much less fre¬ 
quently observed than standard type III bursts and are observed 
at mu ch higher frequencies during flares (e.g. [Melendez et al 


\\999) . This is due to a number of reasons: difficulty of the radio 
waves to escape the low corona due to absorption, an increase in 
the background electron density, and the presence of a positive 
background electron density gradient. Furthermore, it is harder 
to excite Langmuir waves in denser plasma. The growth rate of 
Langmuir waves is inversely proportional (amongst other terms) 
to background density or oc . In a higher density environ¬ 
ment, Langmuir waves are refracted out of resonance with the 
electron beam to higher phase velocities (lower values in wave 
number space or k-space) where a smaller flux of electrons is 
presen t (see e.g. |Kontar||2001 [[Reid & Kont^|2Q13t|Krafft et aT 


[2013[ ). This results in fewer Langmuir waves to convert to radio 
waves and may explain the reason why the reverse drifters stop 
emitting sooner. As an example of this, bi-directional type III 
bursts at high freq uencies (in the metric /decimetric range) have 
been simulated by [Li et aL] ( [2008[[2011[ ). They highlight the en¬ 
hanced difficulty that an electron beam has in producing radio 
emission with a positive drift rate. 


A second factor at play here could be the collisional losses 
experienced by the electron beams. To test this assertion, we 
were able to estimate the column densities encountered by the 
electron beams from f n{r)tcddr, where dr is integrated along 
the path length of the beam and n{r)tcd is from our density 
model. We find that the anti-sunward beams encounter a col¬ 
umn density of 1.4 x 10^^ cm“^, while the sunward beams en- 



Beam velocity (c) 


Fig. 8. (a) Scatter plot of beam displacement (using the TCD 
density model) against burst maximum intensity above the back¬ 
ground, for reverse drift burst only. All bursts together show a 
weak positive correlation coefficient of 0.48, meaning the more 
intense the bursts the further they may travel in space. There is 
a clear separation in intensity between the two types of reverse 
drift bursts. By themselves, the reverse drift bursts starting at 
43 MHz have a correlation coefficient (CC) of 0.60, while those 
starting at 32 MHz have CC of 0.55. The black circles with ar¬ 
rows indicate that these bursts may be slightly longer in space. 
Their corresponding bursts drifted outside the dynamic spectra 
observing range, (b) Maximum intensity vs beam velocity, there 
is almost no correlation here (CC = 0.31), with the only depen¬ 
dency being with the two different sets of bursts, i.e. one set 
(green points - start frequency of 32 MHz) are less intense and 
slower than the other. 


counter 1.1 x 10^^ cm“^. The lower column densities for the out¬ 
ward propagating beams may not be sufficient to significantly 
damp the beams, especially for the beam speeds we observe in 
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Figure 5. However, the average column densities for the down¬ 
ward propagating beams are higher (1.1 x 10^^ cm“^), enough 
to collisionally stop electrons at 3 keV (0.1c) or less ( |Brown[ 
1973t Emslie & Smith[p^84 ). Given the sunward beam speeds 


of 0.14c, it is possible that Coulomb collisional energy losses 
play a role in stopping these beams. However, it is not clear 
whether collisions are the dominant energy loss mechanism that 
causes radio emission to cease. We observe many drift rates with 
higher velocities in Figure 5 and we do not observe a correlation 
between inferred drift rate velocities and inferred beam displace¬ 
ment. We simply highlight collisional losses due to higher col¬ 
umn densities as a possible factor in the shorter displacements 
of the sunward beams. 


5.3. Intensity, velocity, and number density 

We also investigated the dependency of burst intensity on other 
characteristics. A number of theore tical studies have been per- 
formed on radio burst emissivit y ([Melrose |1986t Robinson 


1992t[Robinson et ari|1993||1994||Li et al.l|2Q12t|Ratcliffe et al 


jM(r) : 


Om t^btneVlAv], 
AQm 3/(r) Vb 


( 1 ) 


here the M stands for either fundamental F or harmonic H emis¬ 
sion, AFIm is the solid angle over which the emission is spread, 
Hb is the electron beam number density, Vb is the beam speed, l(r) 
is the distance from emission point to observer, Avb is the width 
of the beam in velocity space, and ^f,m are wave conversion 
efficiency factors which depend on v^. 

However, we find only a weak correlation between burst in¬ 
tensity and velocity (Figure [^). The weak correlation could be 
related to variations in the beam density; a parameter that we 
are unable to infer. If the accelerator energises electron beams 
with a weak correlation between beam density and beam veloc¬ 
ity then we would not expect a tight correlation between burst 
intensity and beam velocity. Moreover, the presence of density 
fiuctuations, known for modulating the type III burst intensity 


(e.g.jMuschietti et ^ [T%5l|Reid & Kontar|[20T0l|Ziebell et al.| 

|2011[|Li et al.[|201^ can increase the scatter between beam ve- 

locity and maximum intensity. 

The acceleration region produces electron beams with more 
intensity when the starting frequency shifts to 43 MHz. Average 
values for the maximum intensity, beam displacement and beam 
velocity all increase (Figure 8) for the reverse drift herring¬ 
bone bursts. Figure 8 suggests that the faster beams that produce 
higher intensity bursts propagate further toward the solar surface 
to reach higher frequency pla sma. This is similar to the results 
of [Cairns & Robins'on| ( |1987| ), in which they showed that bursts 
with greater velocity reached further distances from the accel¬ 
eration site. Interplanetary type Ills display this property with 


2014). The general process of plasma emission involves the for¬ 
mation of a particle beam through some acceleration process; the 
plasma becomes unstable due to the presence of this beam, re¬ 
sulting in the growth of Langmuir waves; these Langmuir waves 
coalesce with ion acoustic waves or secondary Langmuir waves 
to produce electromagnetic emission at the local plasma fre¬ 
quency and its harmonic. To compute the expected emissivity 
of the resulting electromagnetic (radio) waves, stochastic growth 
theory has been employed in a number of cases to successfully 
explain the observed intensities of radio bursts ([Robinson &[ 
[Caims[[l993t[Knock et al.[[2003[[Schmidt & Caims[[2Q12[ ). The 

theory implies that burst emissivity should depend on the beam 
velocity and number density, given by 


Leblanc et al. ( 1995[p^96] ) finding lower stopping fre quencies 
for type III bursts with higher radio flux. Recent work by [Reid &| 
[Kontar[p015[ ) found that beam density and injected spectral in¬ 
dex both affect the stopping frequency of type III bursts. Overall, 
Figure would suggest that beams that are greater in number 
density and/or velocity travel a larger distance in the corona. 


5.4. Accelerator environment change 


Midway through the radio activity, the origin of acceleration ei¬ 
ther drops significantly in altitude from 0.59 R© to 0.47 R© or 
experiences an environment change. This is detectable through 
the increase in start frequency from 32 MHz to 43 MHz of both 
forward and reverse drift bursts. The shift occurs at around 
10:50:30 UT and a change in the characteristics in herringbones 
can be seen. Firstly, the reverse drift bursts are initially morpho¬ 
logically different in the dynamic spectrum, with some bursts 
either showing an increase in intensity over time and others 
showing an increase then decrease. After the shift, the herring¬ 
bones acquire a faster drift rate and larger (more negative) rate 
of change of intensity. Figure [7] shows this clearly, and confirms 
that there is a morphological difference in the herringbones be¬ 
fore and after the shift in altitude. This implies that the presence 
of herringbones in a radio shock signature could have a depen¬ 
dency on the environment in which the shock propagates. This 
is in addition to the theories of a wavy shock front producing 
herringbones ( [Zlobec et al.[ [1993t [Vandas & Karlicky] [2011[ ), 
i.e. herringbone formation may depend as much on coronal en¬ 
vironment as it does on shock characteristics, which has not been 
shown in the past. 


5.5. Alternative acceleration mechanisms? 


Finally, herringbones have traditionally been interpreted as be¬ 
ing formed by a s hock and there is much evidence to sug- 
gest a shock origin ( [Cairns & Robinson| [1987[ [Cane & White 


1989). Various models have been produced to explain the exis¬ 


tence of herringbones, such as a wavy shock front on a CME 
flank ([Zlobec et al.] [1993t [Vandas & Karlicky] [20111 [Schmidt 
[& Cairnsj [2012[) or the termination shocks of a re connection 
outflow jet ( Aurass et al.j [2002[ [Mann et aT) [2009| ). However, 
an explanation of herringbone burstiness is still not complete. 
This bursting has been shown to have a quasi-periodicit y of 2- 
11 seconds, at least for the 2011 September 22 event ( Carley 
[et al.j [2013[ ). In this study we have shown that herringbones 
may have as much to do with the coronal environment as well 
as the shock acceleration mechanism, hence it may be the en¬ 
vironment which produces this bursting. For example, a shock 
which quasi-periodically encounters changes in the coronal en¬ 
vironment which promote the acceleration of particles could pro¬ 
duce such a bursty radio signature. Alternatively, there are other 
mechanisms occurring in the corona which could lead to such 
bursty acc eleration, such as a tearing mode instability in a cur¬ 
rent sheet ( [Kliem et"ari[2000[ ). However, such a mechanism has 
only been attributed to particle acceleration observed at GHz fre¬ 
quencies, much higher frequency (much lower altitude) than the 
bursts we observe here. Further observation and modelling are 
needed to confirm or rule out any of the above mentioned mech¬ 


anisms. 
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6. Conclusion 


Herringbone radio bursts are direct signatures of particle acceler¬ 
ation occurring at coronal shocks. Given the high number of in¬ 
dividual bursts, they provide a good opportunity to perform sta¬ 
tistical analysis on the particle beam properties and/or the prop¬ 
erties of the coronal environment into which the beams propa¬ 
gate. It is because of the high number of individual bursts that an 
automated routine to analyse them quickly and efficiently should 
be developed. 

We have shown here that the Hough transform can provide a 
good automated feature recognition routine to analyse herring¬ 
bones, given that they appear as straight lines in the dynamic 
spectrum. However, the algorithm has some shortcomings, the 
most significant of which is failing to detect the start and stop 
frequencies of the bursts. As a results this step must be done by 
a human user, making the whole process semi-automated. Some 
extra processing steps are needed if all herringbone radio bursts 
are to be analysed for a complete study of their statistics in a 
fully automated way. This is especially pertinent given the latest 
generation of radio telescopes, such as the Low Frequency Array 
(LOFAR; |van Haarlem et al.||2013| ). LOFAR’s frequency range 
of 20-240 MHz, high sensitivity and high frequency and time- 
resolution are ideal for studying the detailed characteristics of 
a variety of fine-structure in s olar radio bursts, such as he rring- 
bones, S-bursts and Type Ills ( [Morosan et al.[ 2014||2015] ) - we 
are set to make great gains in the knowledge of these phenomena 
with current and upcoming radio observational technology. 

We have shown here that a statistical analysis of herringbone 
fine structure can yield useful information on particle beam kine¬ 
matics, the coronal environment that the bursts propagate in, and 
the emission mechanisms that are at play for these kinds of ra¬ 
dio bursts. An analysis of a much larger set of herringbone bursts 
could yield important information in each of these areas, as well 
as produce a definitive theory of herringbone production, which 
does not yet exist. 
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